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Abstract. We report on new attempt to investigate baryon-baryon interactions in lattice QCD. 
From the Bethe-Salpeter (BS) wave function, we have successfully extracted the nucleon-nucleon 
(AW) potentials in quenched QCD simulations, which reproduce qualitative features of modern 
AW potentials. The method has been extended to obtain the tensor potential as well as the central 
potential and also applied to the hyperon-nucleon (YN) interactions, in both quenched and full QCD. 
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INTRODUCTION 

In 1935 Yukawa introduced virtual particles, pions, to explain the nuclear force[l], 
which bounds protons and neutrons inside nuclei. Since then the nucleon-nucleon (A/A/) 
interaction has been extensively investigated at low energies both theoretically and 
experimentally. Fig. 1 shows modern A/A/ potentials, which are characterized by the 
following features [2, 3]. At long distances (r > 2 fm ) there exists weak attraction, which 
is well understood and is dominated by the one pion exchange, while contributions from 
the exchange of multi-pions and/or heavy mesons such as p, CO and o lead to slightly 
stronger attraction at medium distances (1 fm < r < 2 fm). On the other hand, at short 
distances (r < 1 fm), attraction turns into repulsion, and it becomes stronger and stronger 
as r gets smaller, forming the strong repulsive core[4]. The repulsive core is essential 
not only for describing the A/A/ scattering data, but also for the stability and saturation of 
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FIGURE 1. Three examples of the modern AW potential for the l So (spin singlet and 5-wave) state: 
Bonn[7], Reid93[8] and AV18[9]. 



atomic nuclei, for determining the maximum mass of neutron stars, and for igniting Type 
II supernova explosions [5]. Although the origin of the repulsive core must be related to 
the quark-gluon structure of nucleons, it remains one of the most fundamental problems 
in nuclear physics for a long time[6]. It is a great challenge for us to derive the nuclear 
potential including the repulsive core from (lattice) QCD. 

In this talk, we first explain our strategy to extract AW potentials theoretically from 
the first principle using lattice QCD and present the recent result in quenched QCD 
simulations. We then apply the method to various cases including the energy dependence 
of the potential, the quark mass dependence of the potential, the tensor potential, the full 
QCD calculation and the hyperon-nucleon potential. Since we mainly present only the 
results due to the limitation of the length, please see the corresponding references for 
more details. 



STRATEGY TO EXTRACT POTENTIALS IN QCD 

Since a potential is a concept of the non-relativistic quantum mechanics, it is non-trivila 
to define it in QCD. To find a reasonable definition of the AW potentials in QCD, we 
first consider the S-matrix below inelastic threshold of the AW scattering. The unitarity 
leads to 

S = e 2iS (1) 

where the "phase" 8 is a hermitian matrix. We next introduce the equal time Bethe- 
Salpeter (BS) wave function[10, 11], defined by 

<p E {r) = (0|Af(x + r,0)jV(x,0)|2Ar,£> (2) 

where \2N, E) is a two-nucleon eigenstate in QCD with energy E = l^rn^ + k 2 , and the 
N(x) is the gauge-invariant 3-quark operator given by 

N(x) = £ abc q a (x)q h (x)q c (x). (3) 
For large r = |r|, the partial wave / of the BS-wave function behaves as 

<p l E (r) -> Ai — i ± K -LL (4) 

where 8i(k) is the phase of the S-matrix for the partial wave Z[12, 13]. (Although 8i(k) 
is a hermitian matrix in general, we here consider the case that 8i(k) is just a number 
for simplicity.) The above formula says that 8i(k) is the scattering phase shift of the 
scattering wave. In other words, the BS wave function defined above can be interpreted 
as the AW scattering wave. 

Based on the above fact, we have proposed the following strategy to define and extract 
the AW potential in QCD[14, 15]. We define a non-local potential from the BS wave 
function as 

[E-H ](p E (r) = Jd 3 S U(r,s)(p E (s) (5) 
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FIGURE 2. The (effective) central potential for the l S ( 3 Si) state at m K = 529 MeV in quenched QCD. 



where Hq = — V 2 / (2/1^) and /i^ = m^/l is the reduced mass of a two-nucleon system. 
Since the non-local potential is difficult to deal with, we expand it in terms of derivatives 
as U (r,s) = V(r, V)6 3 (r — s). The first few terms are given by[16] 

V(r,V) = Vc(r)+V T (r)S l2 + V L sL-S + {V D (r),V 2 } + ---, (6) 
3 

S12 = - 1 (o 1 -r)(o 2 -r)-(o 1 -o 2 ) (7) 

where S12 is the tensor operator and a ; is the spin operator of the z-th nucleon. The central 
potential Vq and the tensor potential Vt are the leading local terms (without derivatives), 
and thus can be determined from the BS wave function at one energy using eq.(5). Higer 
order terms in the above derivative expansion can be successively determined from BS 
wave functions at different energy. We however expect at low energy that the leading 
order terms, Vc and Vt, give a good approximation of the potential. We can estimate an 
applicable range of energy for the local potential approximation, by calculating physical 
observables such as the scattering phase shift from the local potential and comparing 
them with experimental values. 

The first result for the AW potentials based on the above strategy has been obtained 
in quenched lattice QCD simulations [14], where the lattice spacing a is 0.137 fm, the 
spatial extension L is 4.4 fm and the pion mass m n is 529 MeV. In Fig. 2, the central 
potential for the 1 Sq (spin singlet and L = 0) state and the effective central potential 
for the 3 5i (spin triplet and L = 0) state, obtained at k 2 ~ 0, are plotted as a function 
of r. By comparing Fig. 2 with Fig.l, we see that qualitative features of AW potentials 
are reproduced. Ref.[14] has been selected as one of 21 papers in Nature Research 
Highlights 2007 [17]. 



RECENT DEVELOPMENTS 



In this section we report recent developments of lattice QCD calculations for baryon- 
baryon interactions, based on the method in the previous section. 



Energy dependence 

We first investigate the applicable range of energy for the local potential determined 
at k ~ 0. If terms with derivatives such as Vls ( r )L ■ S or {Vb(r) , V 2 } becomes important, 
the local potential determined at k > is different from the one at k ~ 0. From such 
k dependences of local potentials, in principle, some of the terms with derivatives can 
be determined. In Fig. 3, the local potential for the So state obtained at k ~ 250 MeV 
(red, APBC) is compared with the one at k ~ MeV (blue, PBC) in quenched QCD at 
a = 0.137 fm and m K = 529 MeV[18, 19]. 

As can be seen from the figure, the k dependence of the local potential turns out to 
be very small. This means that the potential obtained at k ~ in Fig. 2 well describes 
physical observables such as the phase shift 6o(k) from k ~ to k ~ 250 MeV, in 
quenched QCD at a = 0. 137 fm and m n = 529 MeV. 



Quark mass dependence 

A quark mass dependence of the NN potential is shown in Fig. 4, where the central 
potentials for the 1 Sq state obtained at k ~ at m n = 380, 529 and 731 MeV are 
compared in quenched QCD at a = 0.137 fm[20, 13]. 

As quark mass decreases, the repulsion at short distance (the repulsive core) get 
stronger while the attraction at intermediate distance ( 0.6 fm ~ 1.2 fm) becomes also 
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FIGURE 3. Comparison of central potentials in the local approximation between the periodic boundary 
condition (PBC, blue) and the anti-periodic boundary condition(APBC, red) for the 1 Sq state at m n = 529 
MeV. 
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FIGURE 4. The central potential for the % state at m % = 380 MeV (red), 529 MeV (green) and 731 
MeV (blue) in quenched QCD at a = 0. 137 fm. 



stronger. 



Tensor potential 

The tensor operator S\i mixes the 3 5i( 7 = 5=1 and L = 0) state with the 3 /J>i 
(7 = 5 = 1 and L = 0) state. Using this property, we can determine the tensor potential 
as follows. For the 7 = 5=1 state, the local potential approximation leads to 

(E-H )<p E (r) = [V c (r) + V T (r)S l2 }<PE(r), (8) 

which, by the projection P to the L = state and the projection Q to the L = 2 state, is 
decomposed into 

/ Pcp E PS 12 (Pe \ v ( Vc\_ (F H J P<Pe\ (Q s 

The above equation can be easily solved as 

(l C ) = (n* n^Y-iE-^f^). (10) 
\V T J \ Q(p E QSn(pE J K >\Q<Pe J v ' 

In Fig. 5, the central potential Vc(r) and the tensor potential Vr(r) for the spin-triplet 
state are plotted[13], together with the effective central potential V^ ff (r) for the 3 5i in 
Fig. 2, which corresponds to (E — Ho)P(p£ / (P<Pe) in the above notation. These potentials 
are calculated in quenched QCD at a = 0. 137 fm and m n = 529 MeV. 

We first notice that the tensor potential Vj has no strong repulsive core, in contrast 
to the central potential. This feature is consistent with the previous phenomenological 
estimate[21]. Although the tensor potential is comparable in the magnitude with the 
central potential, the difference between the central and effective central potentials, 



100 



> 



50 - 



- 



-50 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r- 

\ V T (r) _ 

V c (r) 



0.0 



V c e,f (r) 



m =529 MeV 



f M » i I if f 



0.5 



1.0 
r[fm] 



1.5 



2.0 



FIGURE 5. The central potential(blue) and the tensor potential (red), together with the effective central 
potential (red), for the spin-triplet state at m K = 529 MeV. They are obtained from k ~ states in quenched 
QCD at a = 0.137 fm. 
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FIGURE 6. Quark mass dependence of the tensor potential. 



which is caused by the second order perturbation of Vt, is very small at this quark mass. 
A quark mass dependence of the tensor potential is given in Fig. 6, where the tensor 
potential is plotted at m n = 380,529 and 731 MeV[13]. The tensor potential becomes 
stronger as the quark mass decreases. 



Full QCD calculations 

Results for the AW potentials so far have been calculated in quenched QCD. In 
Ref.[22], preliminary results in full QCD calculations have been reported, based on 
gauge configurations generated by the PACS-CS Collaboration in 2+1 flavor QCD at 
a = 0.09 fm and L = 2.9 fm[23]. Hadron spectra obtained from these configurations are 
shown in Fig.7. Agreements between lattice QCD predictions and experimental values 
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FIGURE 7. Light hadron spectra extrapolated to the physical point using m K , iuk and ihq as input. 
Horizontal bars denote the experimental values. 



3500 
3000 
2500 
2000 
1500 
1000 
500 




150 



100 



50 



-50 



* 



V c (r;'S ) 
Vc'VSi) 



m T =702 MeV 




■ ■ . . . l . , . , l . . . . l_ 



0.0 0.5 1.0 1.5 2.0 2.5 



0.0 



0.5 



1.0 1.5 
r[fm] 

i < 



2.0 



2.5 



FIGURE 8. The (effective) central potential for the So state (red) and the S S\ state (blue) at m K = 702 
MeV in 2+1 flavor QCD at a = 0.09 fm. 



are quiet good. The (effective) central NN potentials on these configurations are given 
in Fig. 8 for the % state (red) and the 3 5i state (blue) at m n = 702 MeV. 

We observe that the repulsive core in both states is much larger in magnitude than 
the corresponding one in quenched QCD in Fig. 2, though the lattice spacing ( 0.09 fm 
vs. 0.137 fm) and the pion mass ( 529 MeV vs. 702 MeV ) are different. A reason 
for the difference of the repulsive core between full and quenched QCD is now under 
investigation. 



Hyperon-Nucleon interactions 

A hyperon is a baryon which contains at least one strange quark. Contrary to the 
case of AW interactions, hyperon-nucleon (YN) or hyperon-hyperon(lT) interactions 
can not be precisely determined, since the scattering experiments are either difficult or 
impossible due to the short life times of hyperons. Our approach therefore may open a 
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FIGURE 9. Left: The central potential for the AN ( l S ) state in 2+1 full QCD as a function of r at 
nii ~ 400 MeV (red) and 700 MeV (green). Right: The central and tensor potentials for the AN ( 3 Si - 3 L>i) 
state in 2+1 full QCD at m K ~ 400 MeV (red and blue) and 700 MeV (green and magenta). 



new possibility to determine them theoretically from QCD. In this direction, potentials 
between a S° (hyperon with strangeness —2 ) and a proton has already been calculated 
in quenched QCD [24]. 

Recently calculations have been extended to the potentials between the A( hyperon 
with strangeness — 1 ) and N in both quenched and full QCD. Fig. 9 shows the AW 
potentials as a function of r obtained from 2+1 flavor QCD calculation[25] based on 
the PACS-CS gauge configurations. The central potential for the x Sq state is given in 
Fig.??, while the central and the tensor potentials for the 3 5i — 3 D\ state are given in 
Fig.?? , with highlighting the short distance (medium to long distance) region in the left 
(right) panel. These figures contain results at m n ~ 400 and 700 MeV. 

As can be seen from both figures, the attractive well of the central potential moves to 
outer region as the quark mass decreases while the depth of these attractive pockets do 
not change so much. The present results show that the tensor force is weaker while the 
spin dependence is stronger than the AW case[22]. As in the AW case, the hight of the 
repulsive core is much larger than the quenched case[25] and it increases as the quark 
mass decreases. 



CONCLUSION 

We present recent results on baron-baryon interactions obtained from lattice QCD simu- 
lations. In our strategy, baryon-baryon(AW, YN and YY) potentials are extracted from the 
BS wave functions. The first result for the AW (effective) central potentials in quenched 
QCD shows good "shape": Qualitative features of the AW potential have been repro- 
duced in quenched QCD, and the energy dependence of the potentials is weak at low 
energy. The method has been successfully extended to the tensor potential and the AW 
potentials in both quenched and full QCD. 

One of the ultimate goal in our approach is to calculate baryon-baryon potentials in 
full QCD at m% =140 MeV. In such calculations one can investigate, for example, a 
relation between the deuteron binding and the tensor force. As other directions of our 
approach, it is important to extract the 3-body force[26] from QCD and to understand 
the origin of the repulsive core theoretically [27]. 
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